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ABSTRACT 

A family of models of thin discs and spheroidal haloes with masses in a linear relationship 
is presented. The models are obtained by considering the gravitational potential as the super- 
position of two independent components, a potential generated by the thin galactic disc and a 
potential generated by the spheroidal halo. The models leads to an expression for the circular 
velocity that can be adjusted very accurately to the observed rotation curves of some specific 
galaxies, in such a way that the models are stable against radial and vertical perturbations. Two 
particular models for galaxies NGC4389 and UG C6969 are o btained by adjusting the circular 
velocity with data taken from the recent paper bv IVerheiien & Sancicii ( f2001l) . The values of 
the halo mass, the disc mass and the total mass for these two galaxies are computed in such 
a way that we obtain a very narrow interval of values for these quantities. Furthermore, the 
values of masses here obtained are in perfect agreement with the expected order of magnitude 
and with the relative order of magnitude between the halo mass and the disc mass. 
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One of the oldest and most important problems in galactic dynam- 
ics is the determination of the mass distribution bas ed on the obser- 
vation s of the circular velocity or rotation curve jPierens & Hurd 
defined as the speed of the stars moving in the galactic plane 
in circular orbits around the center. Now, if we assume a particu- 
lar model for the composition of the galaxy, the fit of that model 
with the rotation curve of a particular galaxy can, in principle, 
completely determine the distribution of mass. So then, the rotation 
curve provides the most direct method to measure the distribution 
pf mass of a galaxy (.Binnev & Tremaine,2008.) . 

Currently, the most accepted description of the composition 
of spiral galaxies is that a significant portion of its mass is concen- 
trated in a thin disc, while the other contributions to the total mass 
of the galaxy comes from a spherical halo of dark matter, a cen- 
tralbulge and, perhaps, a central black hole tBinnev & Tremain3 
I2OO8I) . Now, since all components contribute to the gravitational 
field of the galaxy, obtaining appropriate models that include the 
effects of all parts is a problem of great difficulty. However, the 
contribution of each part is limited to certain distance scales, so in 
a reasonably realistic mo del is not nec essary to include the contri- 
bution of all components ( lFabeil2006l) . 

In particular, the gravitational influence of the central black 
hole is appreciable only within a few parsecs around the center of 
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the galaxy jSchodel et ^l2002h . so it can be completely neglected 
when studying the dynamics in the disc, or in regions outside the 
central bulge, while the bulge mainly dominates the inner region 
of the galaxy to a few kiloparsec. So then, the main contributions 
to the gravitational field o f the galaxy come from the galactic disc 
and the dark matter halo ( lFabeij|2006h . However, it is commonly 
accepted that many aspects of galactic dynamics can be described, 
in a fairly approximate way, usin g models that consider onl y the 
contribution of a thin galactic disc l lBinnev & Tremain3l2008l) . 

Accordingly, the study of the gravitational potential gener- 
ated by an idealized thin disc is a problem of great astrophys- 
ical relevance and so, through the years, different approaches 
have been used to obtain such kind of thin disc models (see 
Isinnev & Tremain3 l l2008l) and references therein). So, once an 
expression for the gravitational potential has been derived, corre- 
sponding expressions for the surface mass density of the disc and 
for the circular velocity of the disc particles can be obtained. Then, 
if the expression for the circular velocity can be adjusted to fit the 
observational data of the rotation curve of a particular galaxy, the 
total mass can be obtained by integrating the corresponding surface 
mass density. 

However, although most of these thin disc models have surface 
densities and rotation curves with remarkable properties, many of 
them mainly represent discs of infinite extension and thus they are 
rather poor flat galaxy models. Therefore, in order to obtain more 
realistic models of flat galaxies, is better to consider methods that 
permit obtaining finite thin disc models. Now, a simple method to 
obtain the gravitational potential, the surface density and the rota- 
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tion cu rve of thin discs of finite radius was developed bv lHunte3 
i ll963h . the simples t exampl e of a disc obtained by this method be- 
ing the well known Kalnajs (1972) disc^ 

In a previous paper ( Gonzalez & Reind l2006h we use the 
Hunter method in order to obtain an infinite family of thin discs 
of finite radius with a well behaved surface mass density. This fam- 
ily of disc models was derived by requiring that the surface density 
behaves as a monotonously decreasing function of the radius, with 
a maximum at the center of the disc and vanishing at the edge. 
Furthermore, the motion of test particles in the gravitational fields 
generated by the first four members of this family was studied in 
iRamos-Car o. Lopez-S uspez &Gonzalez (2008). So, although the 
mass distribution of this family of discs presents a satisfactory be- 
haviour in such a way that they could be considered adequate as flat 
galaxy models, their corresponding rotation curves do not present 
a so good behavior as they do not reproduce the flat region of the 

observed rotation curve. 

On the other hand, in IPedraza. Ramos-Caro &Gonzalezl 

( l2008h a new family of discs was obtained as a superposition of 
members of the previously obtained family, by requiring that the 
surface density be expressed as a well behaved function of the grav- 
itational potential, in such a way that the corresponding distribution 
functions can be easily obtained. Furthermore, besides presenting 
a well-behaved surface density, the models also presented rotation 
curves with a better behavior than the generalized Kalnajs discs. 
However, although these discs are stable against small radial per- 
turbations of disc star orbits, they are unstable to small perturba- 
tions vertical to the disc plane. Then, apart of the stability problems, 
these discs can be considered as quite adequate models in order to 
describe satisfactorily a great variety of galaxies. 

Based on these works, in 

iGonzalez. Plata-Plata & Ramos-Carol boioh were obtained 
some thin disc models in which the circular velocities were 
adjusted to very accurately fit the observed rotation curves of 
four spiral galaxies of the Ursa Major cluster, galaxies NGC3877, 
NGC3917, NGC3949 and NGC4010. These models presented 
well behaved surface densities and the obtained values for the 
corresponding total mass it agrees with the expected order of 
magnitude. However, the models presented a central region with 
strong instability to small vertical perturbations. Now, this result 
was expected as a consequence of the fact that the models only 
consider the thin galactic disc. Therefore, more realistic models 
must be considered that include the non-thin character of the 
galactic disc or the mass contribution of the spheroidal halo. 

In agreement with the above considerations, in this paper we 
will consider a family of models obtained by expressing the grav- 
itational potential as the superposition of a potential generated by 
the thin galactic disc and a potential generated by the spheroidal 
halo, in such a way that the model implies a linear relationship 
between the masses of the thin disc and the spheroidal halo. By ad- 
justing the corresponding expression for the circular velocity to the 
observed data of the rotation curve of some specific galaxies, some 
particular models will be analysed. Then, from the corresponding 
expressions for the disc surface density and the density of the halo, 
estimate values for the total mass of the disc and the total mass of 
the halo will be obtained. The paper is organised as follows. First, 
in Section 2, we present the thin disc plus halo model. Then, in 
Section 3, we obtain the corresponding expressions for particular 
models and, in Section 4, the models are fitted to data of the ob- 
served rotatio n curve of some galaxies o f the Ursa Mayor Cluster, 
as reported in lVerheiien & Sancicil (1200 ih . Finally, in Section 5, we 
discuss the obtained results. 



2 THE THIN DISC PLUS HALO MODEL 

In order to obtain galaxy models consisting of a thin galactic disc 
and a spheroidal halo, we begin considering an axially symmetric 
gravitational potential $ = ^{R, z), where {R, ip, z) are the usual 
cylindrical coordinates. Also, besides the axial symmetry, we sup- 
pose that the potential has symmetry of reflection with respect to 
the plane z = 0, 



(1) 



which implies that the normal derivative of the potential satisfies 
the relation 



-{R, ~z) = ~ — {R,z), 



(2) 



dz ' dz 
in agreement with the attractive character of the gravitational field. 
We also assume that d^/dz does not vanish on the plane 2 = 0, in 
order to have a thin distribution of matter that represents the disc. 

On the other hand, in order to separately describe the thin disc 
and the spheroidal halo, we consider that the gravitational potential 
can be written as the superposition of two independent components 



HR,z) = ^d{R,z) + it>h{R,z), 



(3) 



where ^d{R, z) is the part of the potential generated by the thin 
galactic disc, while <I>h(i?,z) corresponds to the spheroidal halo 
component. The disc component <^d[R,z) must be a solution of 
the Laplace equation everywhere outside the disc. 



0, 



(4) 



while the halo component ^h{R, z) satisfies the Poisson equation 



(5) 



where g{R, z) is the mass density of the halo. 

So, given a potential z) with the previous properties, we 
can easily obtain the circular velocity Vc{R), defined as the velocity 
of the stars moving at the galactic disc in circular orbits around the 
center, through the relationship 



vUR) - R ^ 



(6) 



while the surface mass density of the thin galactic disc is 

given by 



E(i?) 



1 a* 



2nG dz 



(7) 



=0+ 



which it is obtained by using the Gauss law and the reflection sym- 
metry of z). 

Accordingly, in order that the potential of the spheroidal halo 
does not contribute to the disc surface density, we will impose the 
condition 



dz 



= 0. 



(8) 



=0+ 



Furthermore, in order to have a surface density corresponding to 
a finite disclike distribution of matter, we impose boundary condi- 
tions in the form 

dz 



/ 0; R^a, 



(9a) 



=0+ 



dz 



= 0; R> a, 



(9b) 



=0+ 



in such a way that the matter distribution is restricted to the disc 
z — 0, ^ R ^ a, where a is the radius of the disc. 
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In order to properly pose the boundary value problem, we in- 
troduce the oblate spheroidal coordinates, whose symmetry adapts 
in a natural way to the geometry of the model. These coordi- 
nates are related to the u sual cylindrical coordinates by the relation 
( lMorse&Fesbachlll953h 



z = air], 



(10a) 
(10b) 



where ^ < oo and — 1 77 < 1. The disc has the coordinates 
^ = 0, ^ r?^ < 1. On crossing the disc, the r] coordinate changes 
sign but does not change in absolute value. The singular behaviour 
of this coordinate implies that an even function of r; is a continu- 
ous function everywhere but has a discontinuous rj derivative at the 
disc. 

Now, in terms of the oblate spheroidal coordinates, the 
Laplace operator acting over any axially symmetric function 
^{i, Tl) gives 

whereas the boundary condition ^ is equivalent to 



cm 

di 

dri 



= 0, 



5=0 



= 0, 



(11) 



(12a) 



(12b) 



?7 = 



5=0 



0. 



(13a) 



(13b) 



and the boundary conditions l [9at and l |9b| l reduce to 

di 

Moreover, in order that the gravitational potential be continuous 
everywhere, ^{£,,rj) must be an even function of 77, which grants 
also the fulfilment of conditions l ll2bt and l ll3bt . 

Accordingly, by imposing the previous boundary conditions 
over the general solution of the Laplace equation in oblate 
spheroidal coordina tes, we can wri te the gravitational potential of 
the galactic disc as jBatemaJl 19441) 



(14) 



where n is a positive integer, which it defines the model of disc 
considered. Here P2i{ii) are the usual Legendre polynomials and 



Q2i(iC), with Q2i(x) the Legendre functions of 



second kind (see Arfken & Weber 20 0jand, for the Le gendre func- 
tions of imaginary argument, Morse & Fesbactill 19531 . page 1328). 
The coefficients C21 are, in principle, arbitrary constants, but which 
must be specified to obtain any particular model. We will do this af- 
ter, by adjusting the circular velocity of the model with the observed 
data of the rotation curve of some specific galaxies. 

With this expression for the gravitational potential of the disc, 
the surface density is given by 



for the total mass of the disc. Now, it is clear that the surface density 
diverges a t the disc edg e, when rj — 0, unless that we impose the 
condition ( lHuntei|[l963 



J2C2li2l + l)q2t + l{0)P2l{0) 
i=0 

that, after using the identities 

P2n(0) = (-1) 
a2n + l(0) 



0, 



^n( 2n-l)!! 

(2n)!! ' 
(2n)!! 



(17) 

(18a) 
(18b) 



(2n+ 1)!!' 

which are easily obtained from the properties of the Legendre func 
tions, leads to the expression 



Co = ^(-1)' + V2 



(19) 



which gives, through ( I16t . the value of the dies mass Aid in terms 
of the constants C21, with 1^1. 

Now, to properly choose the gravitational potential of the 
spheroidal halo, first we consider the product functions 

^g'm'W, (20) 

where jLambI 19451) 



are the solutions of the differential equation 



d_ 

di 



jU + 1) 



i + e 



(21) 



(22) 



while the associated Legendre functions jArfken & Webed[2005b 

2.^d'=Pj{ri) 



are the solutions of the differential equation 



d_ 

dr) 



(I-'?') 



2,dP^ 



drj 



-i(j + i) 



(23) 



(24) 



where j and k are integers, with j ^ k. 

With the previous expressions, and using the Laplace opera- 
tor i ll It . we can now compute the mass density associated with a 
gravitational potential given by $*^(^,»7). Thus, from the Poisson 
equation ([5j, we obtain 



a2(l + ^2)(l_^2 



(25) 



E(i?) 



1 



2iiaGr) 



with fc 7^ 0, while the density is identically zero when fc = 0. 
Then, as Pj{ri) is finite at the interval — 1 77 ^ 1 and QjH) goes 
to zero when ^ — >■ 00, Q'j{i,ri) properly vanishes at infinity. On 
the other hand, from i23t it is easy to see that rf) diverges at 
77 = ±1 (the z axis) unless that fc ^ 2, while for fc > 2 this function 
vanishes at 77 = ±1. On the other hand, due to the discontinuous 
n character of 77, ^^(C,'?) will be continuous across the disc only 

^6*2,(2/ + I)g2!+i (0)P2i(»7), (15) if — k) is an even number, since then (77) will be an even 
;=o function of 77. 

In agreement with the above considerations, we will take for 
the gravitational potential of the spheroidal a simple combination 
of the product functions $j (5, 77) given by 

^-^C,*?) =-Boo$o(C,'?) + S20$2(C,'?)+B22$i(?,77), (26) 



where, as ^ = 0, 77 = \/ 1 — R?, with R = R/a. Then, by inte- 
grating on the total area of the disc, we find the value 



Co 



(16) 
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where Boo, B20 and ,622 are, in principle, arbitrary constants. 
Then, after using the explicit expressions for q^{£,) and Pf{rj), the 
halo density can be written as 



3B22 
7rGa2 



3 cot ^ ^ — 



e(3C^ + 5) 



(27) 



which is maximum at the disc surface, when ^ = 0, and then fastly 
decreases being constant at the oblate spheroids defined by ^ = 
cte. By integrating over all the space, we obtain the expression 



B22 



MhG 
16a ' 



(28) 



where A^h is the total mass of the halo. Finally, from the condition 
J12at . we obtain 



Boo = 
B20 = 



MhG 

a 

MhG 
2a ' 



(29a) 
(29b) 



and so all the constants in l l26t are expressed in terms of the halo 
mass Mh- 

On the other hand, if we restrict to particles moving in the thin 
disc, the circular velocity l|6ll is written in terms of the spheroidal 
coordinates as 



(rf - 1) a$ 



dri 



(30) 



which, by using (O, dl4b . l l26t and the properties of the Legendre 
functions, reduces to 



ri •'^ — ' 

' /—I 



where 



and 



C2 = 92(0) 



C2 + 



MhG 



4a 



(31) 



(32) 



where the A21 constants are related with the previous constants C21 , 
for / 7^ 0, through the relation 



Co, 



4; + 1 
4Z(2Z + 1; 



where 



/: 



77(1 - ■<i)-f P'2i{j])dr), 



(36) 



(37) 



which is obtained by equaling expressions | |3U and l l35t and by 
using t he orthogonality prope rties of the associated Legendre func- 
tions tArfken & Webeij|2005h . 

Then, if the constants A21 are determined by a fitting of the 
observational data of the corresponding rotation curve, the corre- 
sponding values of the coefficients C21 can be determined by means 
of relation l l36t , obtaining then a particular case of l l34t correspond- 
ing to a specific galaxy model, which can be written in terms of the 
constants A21 as 



MiG ^ MhG (- 
;=i fc=i 



{41 + l)A2kIki . 



a 



4a 



41(21 + I)q2ii0) 



(38) 



however, this relation does not determine completely the values of 
Md and Mh, only gives a linear relationship between them. So, in 
order to restrict the allowed values of these masses, it is needed to 
analyse the behaviour of some other quantities characterizing the 
kinematics of the model. 

Besides the circular velocity, two other important quantities 
concerning the kinematics of the models are the epicycle or ra- 
dial frequency, k^{R), and the vertical frequency, u^{R), which 
describe the stability against radial and vertical perturbations of 
particles in quasi-circular orbits (Binnev & Tremaine 2008, Section 
2.2.3, page 164). These two quantities, which must be positive in 
order to have stable circular orbits, are defined as 



(39) 



G21 ^ q2i{0)C2i, (33) 
for / ^ 2. Then, by using l ll9t . l l32t and l l33t . it is easy to see that 

{-iy+'C2i 



MdG ^ MhG ^ ^ 



a 



4a 



92! (0) 



(34) 



and thus the model implies a linear relationship between Md and 
A^h, in which the independent term is determined by the constants 
C21, with 1^1. Now, it is clear that the above relationship makes 
sense only if the right hand side it is positive, which should be 
checked for every set of constants C21 corresponding to any partic- 
ular model. 



where 



(40) 



(41) 



is the effective potential and £ = Rvc is the specific axial angular 
momentum. Then, by using expression ([6j for the circular velocity, 
we can write the above expressions as 



2, , 1 , 2«c 
^^^^^RdR+-R^ 



(42) 



3 OBTAINING OF PARTICULAR MODELS 

In order to obtain particular models, we must to specify the con- 
stants C21 of the general model. So, we will adjust these constants 
in such a way that the circular velocity (R) it fits the data of the 
rotation curve of some particular galaxy. Now, as expression l l31t 
for the circular velocity only involves derivatives of the Legendre 
polynomials of even order, it is equivalent to 

n 
! = 1 



V>|„., 4^, (43. 

where we also used the expression for the Laplace operator in cylin- 
drical coordinates. 

Now, by using d35t . the epicycle frequency can be cast as 

n 

= ^2(« + l)yl2ii?''"', (44) 
(=1 

where k — an. So, as the epicycle frequency is completely deter- 
mined by the set of constants A21, which are fixed by the numerical 
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fit of the rotation curve data, there are not free parameters that can 
be adjusted by requiring radial stabiUty. However, relation l l44t can 
be used as a consistency criterion to the numerical fit of the cir- 
cular velocity with the observed data of the rotation curve. So, if 
the epicycle frequency it is not positive everywhere at the interval 
^ i? ^ 1 for a given set of constants, it is possible to make a new 
numerical fit by considering a different model, with another value 
of n, in order to ensure the positiveness of the epicycle frequency. 
Accordingly, we seek for a balance between the best numerical fit 
and the numerical fit that agrees with the positiveness of the epicy- 
cle frequency. 

Also, by using the Poisson equation J5](, the expression l l27b 
for the halo density and the expression l l35t for the circular velocity, 
the vertical frequency can be written as 

^2,^. 'd-KGMh r ru\ rA'^\ 

^ (R) = 3 /i(-R)' (45) 



8a 



where v = au and 



/i(7?) = ^2M2ii?' 



(46) 



So, as must be positive everywhere at the interval ^ 1 
in order to have vertically stable models, expression l l45t defines a 
lower bound for the halo mass. 



Mh^ M,,= ^ |max/i(i?)| 



(47) 



that combined with ( 138b gives also an upper bound for the disc mass 

We also need to consider the behavior of the surface mass den- 
sity, as given by expression ( llSb . that can be written in terms of the 
constants A21 as 



where, for the model with n = 4, 



(48) 



+ 



2A2 + - 
16^4 



8A4 16^16 128^8 
9 ^ 25 245 



^ 64^ ^ 768^8 



75 

128^6 ^ 1024^8 



75 



1225 



1225 



"2048^8 



1225 



R\ (49) 



and for any other value of n the explicit expression of f2{R) can 
be easily obtained from l ll5t . i32l . i33\ and l l36t . Accordingly, in 
order that the surface mass density be positive everywhere on the 
disc, expression l l48t defines an upper bound for the halo mass, 

Mh^M+ = ^ {mm f2{R)], (50) 

that combined with i38t gives also a lower bound for the disc mass 
Md- However, in order that the models make sense, we must have 
that 



(51) 



which should be checked for every set of constants corresponding 
to any particular model. 



4 FITTING OF DATA TO THE MODELS 

In order to check the applicability of the previous model to describe 
real galaxies, we analyse a set data taken from the recent paper by 



Table 1. Galactic radius a (kpc) and constants A21 (km s ^) 





a 


A2 


A4 


Ae 


As 


NGC4389 
UGC6969 


5.5 
4.6 


29495.4 
15094.2 


-53607.4 
-38417.9 


60170.1 
56183.2 


-23962.6 
-26619.2 



Table 2. Halo mass A^^, disc mass. Aid and total mass Ai (kg). 





Mh (10^4 kg) 


Md (10^3 kg) 


M (1034 kg) 


NGC4389 
UGC6969 


4.24 - 4.52 
1.81 - 1.83 


3.01 - 3.70 
1.42 - 1.47 


4.62 - 4.82 
1.96 - 1.98 


Table 3. Halo 


mass A4h, disc mass. Add and total 


mass M (Mq). 






Aid (Mq) 


M (Mq) 


NGC4389 
UGC6969 


21345 - 22738 
9116 - 9219 


1514 - 1862 
712 - 738 


23207 - 24252 
9854 - 9931 



IVerheiien & Sancicil tOOlh . in which are reported the results of an 
extensive 21 cm-line synthesis imaging survey of 43 galaxies in the 
nearby Ursa Major cluster using the Westerbork Synthesis Radio 
Telescope. For each rotation curve data, we take the value of a, the 
radius of the galactic disc, as given by the last tabulated value of the 
radius. Accordingly, we are assuming that the radius of the galaxies 
is defined by the farthest observed data. Now, although this assump- 
tion about the galactic radius do not ag rees with the accepted stan - 
dard about the edge of the stellar disc dBinnev & Merrifiel3ll998h . 
we will make it since we are assuming that all the stars moving in 
circular orbits at the galactic plane are inside the disc and that there 
are no stars moving outside the disc. 

So, by taking the radius normalized in units of a, we make a 
non-linear least square fit of t he data with the general rela tion i35l 
for all the galaxies reported bv lVerheiien & Sancicil ( l200lh . seeking 
for a balance between the best numerical fit and the fit that agrees 
with the positiveness of the epicycle frequency by using the consis- 
tency criterion ( 144b over the constants A2i. By selecting between 
all the galaxies only those that fulfill with the consistency criterion, 
then we proceed to determine the lower and upper bounds over the 
halo mass by using expressions ( 147b and dSOb . Finally, we check 
that this two values of the halo mass are in agreement with the con- 
dition ( 151b in order that the models make sense. However, when 
we check the consistency of the adjust, we found that only the fit 
of the data for the galaxies NGC4389 and UGC6969 it agrees with 
this conditions, whereas that for all the other galaxies we found that 

In Table [T] we present the values of the galactic radius a and 
the constants A21 obtained by the numerical adjust with the rotation 
curve data for galaxies NGC4389 and UGC6969 by taking a model 
with n = 4. With this values for the constants, we obtain for the 
galaxy NGC4389 the relationship 



Mh + AMd = 5.7255 x 10''* kg 
and for the galaxy UGC6969 the relationship 

Mh + 47Wd = 2.3995 x 10^* kg 



(52) 



(53) 



We also compute the minimum and maximum values of the halo 
mass Mh, the disc mass Md and the total mass M, defined as 



M = Mh + Md 



(54) 
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NGC4389: Circular Velocity 




UGC6969: Circular Velocity 




0.2 0.4 0.6 0.8 1 



Figure 1. Plots of the circular velocity Vc in km/s, as functions of the 
dimensionless radial coordinate R = R/a, for the galaxies NGC4389 and 
UGC6969. 



So, in Table |2l we present the values for the halo mass, the disc 
mass and the total mass in kg, whereas that in Table [3] we present 
the corresponding values in solar mass units [Mq). 

In Fig. [T] we show the adjusted rotation curve for these two 
galaxies. The upper plot corresponds to the model for galaxy 
NGC4389, while the lower plot to the model for galaxy UGC6969. 
The points with err or bars are the observations as reported in 
IVerheiien & SancicH (^001), while the solid line are the rotation 
curves determined from l l35t with n = 4, and the A21 given by 
the numerical fit. As we can see, with the two galaxies we get a 
fairly accurate numerical adjustment to the rotation curve. We also 
plot, in Fig. 12] the epicycle frequency for the two galaxies, and we 
can see that this quantity is always positive in both of them. 

In Fig.|3]we present the vertical frequency for the two models, 
the model for galaxy NGC4389 in the upper plot and the model for 
galaxy UGC6969 in the lower plot. In both of the models we plot 
the vertical frequency for the minimum value M'^ of the halo mass 
(upper curve of each plot) and the maximum value A^^ of the halo 
mass (lower curve of each plot). As we can see, for the two galaxies 
we obtain a positive vertical frequency for all the disc extension, so 
the models are stable against vertical perturbations. Now then, for 
galaxy UGC6969 the two curves are almost coincident given that 
the interval between the two mass values is very narrow. 

Finally, in Fig. [4] we present the corresponding plots of the 
surface mass density for the two galaxies, again with the model 
for galaxy NGC4389 in the upper plot and the model for galaxy 
UGC6969 in the lower plot, plotting the surface mass density for 
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Figure 2. Plots of the epicycle frequency X 10 ^ in (km/s)'^, as func- 
tions of the dimensionless radial coordinate R = R/a, for the galaxies 
NGC4389 and UGC6969. 

the same two values of the halo mass as in the plots of the verti- 
cal frequency. For galaxy NGC4389, the density has a maximum 
at the disc centre and then decreases as 7? increases, vanishing at 
the disc edge. On the other hand, for galaxy UGC6969 the density 
behaves in a similar way, having a maximum at the disc centre and 
then decreasing, presenting then a small region where the density 
remains almost constant and then decreasing again to vanish at the 
disc edge. Furthermore, as with the vertical frequency, for galaxy 
UGC6969 the two curves are almost coincident again as a conse- 
quence of the narrow interval between the two values of the halo 
mass. 



5 DISCUSSION 

We presented a family of galactic models characterised by a lin- 
ear relationship between the halo mass and the disc mass. The 
models are stable against radial and vertical perturbations and 
their circular velocities can be adjusted very accurately to the ob- 
served rotation curves of some specific galaxies. The here pre- 
sented models are a generalisation of the models presented in 
'Gonzalez. Plata- Plata & Ramos-Carol ( l201(]h . where we only con- 
sidered models with a thin galactic disc. The generalisation was 
obtained by adding to the gravitational potential of the thin disc the 
gravitational potential corresponding to a spheroidal halo, in such a 
way that we solve the problem of vertical unstability presented by 
the previous models. 
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Figure 3. Plots of the vertical frequency X 10 ^ in (km/s)^, as func- 
tions of the dimensionless radial coordinate R = R/a, for the galaxies 
NGC4389 and UGC6969. 



Two particular models were obtained by a numerical fit of the 
general expression l l35t for the circular velocity with the observed 
data of the rotation curve of galaxies NGC4389 and UGC6969. For 
these two galaxies we obtained a fairly accurate numerical adjust- 
ment with the rotation curve and, from the constants A21 obtained 
with the numerical fit, we compute the values of the halo mass, the 
disc mass and the total mass for these two galaxies in such a way 
that we obtain a very narrow interval of values for these quanti- 
ties. Furthermore, the values of masses here obtained are in perfect 
agreement with the expected order of magnitude and with the rela- 
tive order of magnitude between the halo mass and the disc mass. 

However, although we tested the applic ability of the present 
model with all the galaxies reported by IVerheiien & Sancica 
( l200lh . consistent models were obtained only for the two galax- 
ies NGC4389 and UGC6969, whereas for all the other galax- 
ies were obtained models with values of the halo mass such that 
Mt < Mn- Now, this negative result can be considered as due to 
the simple halo model that we have taken here. Indeed, as we can 
see from expressions l l25t and i26l . only one term of the gravita- 
tional potential of the halo contributes to their density, what leaves 
only one free constant to be determined in order to fit the model to 
the imposed consistency conditions. This constant is precisely the 
mass of the halo, Aih, which is determined by requiring the posi- 
tiveness of the vertical frequency and the surface mass density. On 
the other hand, by considering additional terms in expression ( 126b 
for the halo potential, we will have new free parameters that per- 
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Figure 4. Plots of the surface mass density E X 10 ^ in (kg/m^), as 
functions of the dimensionless radial coordinate R = R/a, for the galaxies 
NGC4389 and UGC6969. 



haps allow to adjust the model to properly describe the behavior of 
other galaxies besides the two considered here. 

In agreement with the above considerations, we can consider 
the simple set of models here presented as a fairly good approxi- 
mation to obtaining quite realistic models of galaxies. In particular, 
we believe that the values of mass obtained for the two galaxies 
here studied may be taken as a very accurate estimate of the upper 
and lower bounds for the mass of the galactic disc and for the mass 
of the spheroidal halo in these two galaxies. Accordingly, we are 
now working on a more involved model, obtained by including ad- 
ditional terms in expression i26\ for the halo potential, in order to 
get some particular models that can be properly adjusted with the 
observed data of the rotation curve of some other galaxies besides 
the two here considered. 



ACKNOWLEDGMENTS 

GAG and JIR want to thank Vicerrectorfa de Investigacion y Ex- 
tension, Universidad Industrial de Santander, by the financial sup- 
port. JI thanks support from the Spanish Ministry of Science and 
Innovation (grant FIS2010-15492). Also GAG is grateful for the 
warm hospitality of Departamento de Ffsica Teorica, Universidad 
del Pais Vasco, where a main part of this work was made. 



© 0000 RAS, MNRAS 000, 000-000 



8 G. A. Gonzalez, J. Ibdnez and J. I. Reina 



REFERENCES 

Arfken G., Weber H., 2005, Mathematical Methods for Physicists. 
6th Ed. Academic Press. 

Bateman H., 1944, Partial Differential Equations. Dover. 

Birmey J., Merrifield M., 1998, Galactic Astronomy. Princeton 
University Press. 

Binney J., Tremaine S., 2008, Galactic Dynamics. 2nd ed. Prince- 
ton University Press. 

Faber T., 2006, Galactic halos and gravastars: static spherically 
symmetric spacetimes in modern general relativity and astro- 
physics, M. Sc. Thesis in Applied Mathematics, Victoria Uni- 
versity of Wellington. 

Gonzalez G., Plata-Plata S., Ramos-Caro J., 2010, MNRAS, 404, 
468. 

Gonzalez G. A., Reina J. L, 2006, MNRAS 371, 1873. 

Hunter C, 1963, MNRAS, 126, 299 

Kalnajs A. J., 1972, ApJ, 175, 63. 

Lamb H., 1945, Hydrodynamics. 6th Ed. Dover. 

Morse P. M., Fesbach H., 1953, Methods of Theoretical Physics. 

Mc Graw Hill, New York. 
Pedraza J. R, Ramos-Caro J., Gonzalez G. A., 2008, MNRAS 390, 

1587. 

Pierens A., Hure J-M. 2004, ApJ, 605, 179. 
Ramos-Caro J., Lopez-Suspez R, Gonzalez G. A., 2008, MNRAS 
386, 440. 

Schodel R. et al., 2002, Nature, 419, 694. 

Verheinjen M. A. W, Sancici R., 2001, A&A, 370, 765. 



